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Protein-surface interactions are ubiquitous in biological processes and bioengineering, yet are not fully under¬ 
stood. In biosensors, a key factor determining the sensitivity and thus the performance of the device is the 
orientation of the ligand molecules on the bioactive device surface. Adsorption studies thus seek to determine 
how orientation can be influenced by surface preparation, varying surface charge and ambient salt concen¬ 
tration. In this work, protein orientation near charged nanosurfaces is obtained under electrostatic effects 
using the Poisson-Boltzmann equation, in an implicit-solvent model. Sampling the free energy for protein 
GB1D4' at a range of tilt and rotation angles with respect to the charged surface, we calculated the prob¬ 
ability of the protein orientations and observed a dipolar behavior. This result is consistent with published 
experimental studies and combined Monte Carlo and molecular dynamics simulations using this small protein, 
validating our method. More relevant to biosensor technology, antibodies such as immunoglobulin G are still 
a formidable challenge to molecular simulation, due to their large size. With the Poisson-Boltzmann model, 
we obtained the probability distribution of orientations for the iso-type IgG2a at varying surface charge and 
salt concentration. This iso-type was not found to have a preferred orientation in previous studies, unlike the 
iso-type IgGl whose larger dipole moment was assumed to make it easier to control. Our results show that the 
preferred orientation of IgG2a can be favorable for biosensing with positive charge on the surface of 0.05C/m^ 
or higher and 37mM salt concentration. The results also show that local interactions dominate over dipole 
moment for this protein. Improving immunoassay sensitivity may thus be assisted by numerical studies using 
our method (and open-source code), guiding changes to fabrication protocols or protein engineering of ligand 
molecules to obtain more favorable orientations. 

I. INTRODUCTION 

Protein adsorption plays a key role in many biotech¬ 
nological applications, particularly biomaterials and tis¬ 
sue engineering, biomedical implants and biosensors. 

Yet, despite their importance, the specific mecha¬ 
nisms governing protein-surface interactions are not fully 
understood.ll^^ 

In the field of biosensors, protein adsorption needs to 
be engineered to obtain a successful device. Biosensors 
detect specific molecules using a nanoscale sensing ele¬ 
ment, such as a metallic nanoparticle or nanowire cov¬ 
ered with a bioactive coating. The prevailing method to 
modify a sensor surface is via a self-assembled monolayer 
(sam) of a small charged group, with ligand molecules 
layered on top to achieve the desired function. Anti¬ 
bodies are a common choice for the ligand molecules, 
although the newest devices use single-domain or single¬ 
chain fragment molecules.l^ Sensing occurs when a tar¬ 
get biomolecule binds to the ligand molecule, changing 
some physical parameter on the sensor, such as current 
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in nanowires or plasmon resonance frequency in metallic 
nanoparticles. 

One of the factors crucially affecting biosensor perfor¬ 
mance is the orientation of ligand moleculesThese 
have specific binding sites, which need to be accessi¬ 
ble to the target molecule for the biosensor to function 
well. Probing protein orientation is thus one key goal 
of adsorption studies. The aim of this study is to as¬ 
certain how orientation can be influenced by fabrication 
conditions regarding surface preparation, such as surface 
charge and ambient salt content. We consider in particu¬ 
lar the antibody immunoglobulin G near a solid surface at 
different charge concentrations and ionic strengths. Us¬ 
ing a smaller molecule (protein GBID4'), we could first 
confirm agreement of our results with published works 
reporting experiment^ and simulations with a combined 
Monte Carlo and molecular dynamics method.l^ These 
previous works, among others, also concluded that elec¬ 
trostatic interactions are the dominant effect in the orien¬ 
tation of adsorbed proteins. In the case of immunoglob¬ 
ulin G (igG), the protein is relevant for biosensor appli¬ 
cations, but its large size would make all-atom molecular 
simulations quite cumbersome and expensive. For this 
reason, other researchers have studied adsorption of IgG 
using a coarse-grained model that considers each residue 
as a sphere (united-residue model) ,1^ finding that elec¬ 
trostatics dominates the orientation for higher surface 
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charges and that a positive charge can result in the de¬ 
sired “tail-on” placement for the IgGl iso-type, at low 
enough salt concentration. Here, we investigate the pre¬ 
ferred orientations for the IgG2a variant, which other re¬ 
searchers found hard to control. In addition to obtaining 
the preferred orientation at different conditions of charge 
and ionic strength, we also take a detailed look at the 
probability distribution in the parameter space. 

Our model for protein-surface interactions uses the 
Poisson-Boltzmann equations in their integral formu¬ 
lation, representing the protein geometry as a dielec¬ 
tric interface in an implicit solvent. We recently veri¬ 
fied the model against an analytical solution valid for 
spherical geometries and studied its numerical conver¬ 
gence in detail.^ Previous studies on protein-surface in¬ 
teraction using the Poisson-Boltzmann equation showed 
that such a model is adequate as long as c onforma- 
tional changes in the protein are slight ,^2^ and also 
that van der Waals effects can be neglected for realistic 
molecular geometries.l^^ Conformational changes of the 
biomolecule can be ignored in this case because bind¬ 
ing sites need to remain nearly unmodified during the 
biosensor fabrication process.A continuum framework 
has been used in the past to study protein orientation,^^ 
but it included ions explicitly. Other studies have used 
a coarse-grained mod el of the molecule, represented as 
a set of sp heres ,^^^^ assigned effective charges at the 
residue le vel or made approximations to account for 
pH effect 

The sensor element (functionalized with the SAm) is 
represented in our model as a charged surface that in¬ 
teracts electrostatically with the biomolecule. A param¬ 
eter sweep of the protein’s rotation and tilt angles with 
respect to the solid surface provides energy landscapes, 
where the probability of finding the system in a given 
micro-state depends on the total free energy. The con¬ 
tinuum approach can thus provide insights to the condi¬ 
tions (surface charge and salt concentration) conducive 
to a favorable orientation of large proteins, too large for 
all-atom molecular simulation with today’s computing 
power. It can also represent solid surfaces of any geome¬ 
try, and we expect that it may in future assist in the de¬ 
sign of better ligand-molecule immobilization techniques 
for high-sensitivity biosensors. 


II. IMPLICIT-SOLVENT MODEL FOR PROTEINS 
NEAR CHARGED SURFACES 


The implicit-solvent model describes a molecular sys¬ 
tem as a set of continuum dielectric regions, and com¬ 
putes the mean-field potential using electrostatics. For 
the case where a protein is dissolved in a solvent, we re¬ 
quire two of such regions: inside and outside the protein, 
interfaced by the solvent-excluded surface (SES). The SES 


determines the closest a water molecule can get to the 
protein, and we generate it by rolling a spherical probe 
of the size of a water molecule around the protein. The 
dielectric constant inside the protein is low (e = 2 to 4) 
and there are point charges placed at the atomic loca¬ 
tions. The solvent region has the dielectric constant of 
water e ~ 80, and we need to account for the presence of 
salt. This model results in a system of partial differential 
equations where the Poisson equation describes the elec¬ 
trostatic potential inside the protein, and the linearized 
Poisson-Boltzmann equation applies outside the protein. 
On the SES, appropriate interface conditions ensure the 
continuity of the potential and electric displacement. 
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FIG. 1: Sketch of a molecule interacting with a surface: 
fti is the protein, the solvent region, Fi is the SES 
and r 2 a surface with imposed charge. 

In this work, we use an extension of the implicit-solvent 
model to consider the effect of charged surfaces. Such 
is the case of the setup sketched by Figure which is 
described mathematically by the following equations: 


V^0i(r) = — —S{t, Tk) in solute (f2i). 


V^(/)2(r) = K^4>2{r) 
= 02 

901 _ 902 

9n 9n 

902 

-^ 2 ^— = CFq 

on 


in solvent {^ 2 )^ 

on interface Fi, and 
on surface r 2 (1) 


where 0^ is the electrostatic potential in region which 
has a permittivity e^, and ao is a prescribed charge on 
the surface. The surface r 2 could correspond to a device 
such as a biosensor. 

Boundary integral formulation — We apply Green’s 
second identity to the system of partial-differential equa¬ 
tions in Q, and evaluate the resulting equations on Fi 
and r 2 to obtain the following system of integral equa¬ 
tions: 
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The function (t)i^Tj = 0i(i*rj ) is the electrostatic potential 
at a point that approaches the surface Vj from the region 
and K and V are known as the single- and double¬ 
layer potentials, correspondingly: 




y^k 

^L/Y 



£ ^ [G'L/Y(rrfc,rr,)] 4>i,r^ dr, 
^ ^^*.ri^L/y(’'rfc,rr,)dr. (3) 


In Equation Q, Gl and Gy are the free-space 
Green’s functions of the Poisson and linearized Poisson- 
Boltzmann equations, respectively. The single-layer po¬ 
tential of a distribution i/; on a surface T evaluated at r, 
E^(0r), can be interpreted as the potential on r due to 
a charge distribution ip on T. Similarly, K^{'ip-p) can be 
seen as the potential induced by a double layer of charges 
(0) with opposite sign at T. 

Rearranging terms, we write Equation © in matrix 
form, as follows: 
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The boundary-integral formulation is not limited to 
represent the protein with a single surface, but can ac¬ 
count for solvent-filled cavities inside the protein region 
and Stern layer s.^^ In those cases, more than one surface 
is required to appropriately represent the protein. Our 
implementation follows the guidelines from Altman and 
CO- worker^ to deal with multiple surfaces. 

The boundary-integral formulation of the implicit- 
solvent model is a popular alternative to compute sol¬ 
vation energies of proteinsbut the effect of charged 
surfaces has rarely been considered. The only work that 
we know of that does include these effects is limited to 
plane, infinite surfaces .1^ 


^0i,ri I = — r Ti, 

.9n J ei ^47r|rr,-rfcl 

+ =0 onTi, 

- Kp(,/.2.rJ + Vy = 0 on r2. 


( 2 ) 


III. METHODS 


A. Discretization and implementation details 


We solve the system in Q numerically using a bound¬ 
ary element method (bem). To represent the SES, we use 
flat triangular panels where the potential (0) and its nor¬ 
mal derivative (dcp/dn) are constant, and then collocate 
the discretized equation on the center of each panel. This 
transforms the integral operators in the matrix equation 
into block matrices of size Np x Ap, where Np is the 
number of panels. Each entry of the block matrix is an 
integral over one panel (Tj), evaluated on the center of 
panel T^: 


VL,ij=[ Gi(rr„rr,)dr,-. (5) 


We classify the integrals in Equation in three 
groups, depending on the distance d between the panel 
and the collocation point. When the collocation point is 
inside the panel being integrated, we get a singular inte¬ 
gral that we solve with a semi-analytical approaclP^ plac¬ 
ing Gauss nodes on the sides of the triangle. We call near - 
singular integrals those where d < 2L {L = V2 • Area). 
Eor near-singular integrals, we use a high-order Gauss 
quadrature rule with 19 or more nodes. Einally, when 
the panel and collocation points are further than 2L from 
each other, we only need 1, 3, 4 or 7 Gauss nodes per el¬ 
ement to get good accuracy. 

To solve the resulting linear system, we use a gen¬ 
eral minimal residual method (gmres). The most time 
consuming part of the GMRES solver is a matrix-vector 
multiplication—in principle, an 0{N^) operation—done 
within every iteration of the solver. But by using 
a treecode algorithm, we perform this operation in 
0{N log N) time.l^ More details on our implementation 
of the BEM can be found in our earlier work,l^ and a 
companion paper 
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B. Energy calculation 

We can decompose the total free energy into Coulom- 
bic, surface, and solvation energy: 

-^Total ~ -^Coulomb “1“ surf “1“ F golv (b) 

Coulomhic energy — The Coulombic energy arises 
simply from the Coulomb interactions of all point 
charges. We compute it by 

i^couiomb = 2 IZ IZ 47r|ri-r,-| 

j i ' ^ ' 

Solvation free energy — The solvation energy is the 
energy contribution of the protein’s surroundings: sol¬ 
vent polarization, charged surfaces, and other proteins. 
We compute it as 
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FIG. 2: Setup of the problem for our 
orientation-sampling studies. 


-^solv — o / ^ (0total ^Coulomb) (8) 

^ Jn 

— ^ ^ Qkjf^total ^Coulomb)(l*/c) 5 (9) 

k=0 

where p is the charge distribution, consisting of point 
charges (which transforms the integral into a sum), and 

0reac ~ 0total ^Coulomb Is 

</>reac.r. = (<Al,r J + Vl'^ ( 10 ) 

Surfaee free energy — We use the description of free 
energy of a surface with prescribed char ge (like r 2 in 
Figure Q from Chan and co-workersThey describe 
the free energy on a surface as 

F,^,i=tJ^Ge4dT, (11) 

where (j) = GcO-q. 


C. Orientation sampling of a protein near a charged 
surface 


where A is the ensemble of all micro-states, is the 
Boltzmann constant and T the temperature. To obtain 
a probability distribution, we assume that electrostatic 
effects are dominant and use Equation (12), sampling 
^totai for different orientations. We define the orientation 
using the angle between the dipole moment and surface 
normal vectors as a reference (tilt angle), varying from 0° 
to 180°. For each tilt angle, we rotate the protein about 
the dipole moment vector by 360° to examine all possible 
orientations. This process is sketched in Figure 

In this case, micro-states are defined by the tilt (o^tiit) 
and rotational ((arot) angles, and we rewrite the integral 
in the numerator of Equation (12) as: 


f ( ^total A T ^ f f ( ^total \ . 

A “ 7 y ““ 

(13) 

where micro-state A is a range of angles a^ot and amt. la 
biosensors, the ligand is adsorbed on the surface (usually 
covalently), hence we are interested on the orientation of 
the molecule very close to the surface, and don’t consider 
configurations away from it. 


D. Structure preparation 


We are aiming to investigate the orientation of pro¬ 
teins near self-assembled monolayers (sam), specifically 
for biosensing applications. In the framework of the 
implicit-solvent model, we can represent the SAM as a sur¬ 
face charge density, and use Equation @ to compute the 
electrostatic potential. According to the Boltzmann dis¬ 
tribution, the probability of finding the system in micro¬ 
state A depends on the total free energy, Ftotab follows: 


P(A) 


/Aexp(-fg^) dA 
Xvexp(-f^)dA 


( 12 ) 


To assess the adequacy of the implicit-solvent model 
for investigating protein-surface interactions, we stud¬ 
ied the orientation of protein GB1D4' mutant near a 
charged surface, since there are results available in the 
literature that we could compare to: both experimental 
observation^ and simulations using a combined Monte 
Garlo and molecular dynamics approach.^^Eigurej^ shows 
the structure of protein GB1D4' (PDB code IPGB), to 
which we applied mutations E19Q, D22N, D46N and D47N 
to obtain the D4' mutant, using the SwissPdb Viewer 
software.l^ We then studied the orientation of antibody 
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FIG. 3: Structure of protein GB1D4' (PDB code IPGB). 


ter, with 32 GPU nodes featuring dual Intel Xeon E2620 
2.0GHz 6-core processors with dual NVIDIA K20 and 
128 GB of memory. All runs were serial: single-CPU 
and single-GPU. We obtained the van der Waals radii 
and charge distribution using pdb2pqil^ with an amber 
forcefield, and generated the meshes using the free msms 
software.^ In these tests, we did not consider a Stern 
layer for either the protein or the charged surface, nor 
the presence of solvent-filled cavities inside the protein. 



FIG. 4: Structure of immunoglobulin G (PDB code 
lIGT). 


immunoglobulin G iso-type IgG2a (PDB code lIGT), a 
widely used protein in biosensors, whose structure is 
shown in Figure This is a more interesting case from 
the point of view of our application, yet it is more diffi¬ 
cult to study with molecular simulation due to its size. In 
both cases, the vector orientation of the dipole moment 
(used as reference for the tilt and orientation angles) was 
obtained using the location of the point charges at the 
locations of the atoms. 


IV. RESULTS 

The results detailed in this section were obtained using 
an extension of the open-source code PyGBe,!^ account¬ 
ing for the presence of surfaces with imposed charge or 
potential.^ We ran the calculations for protein GB1D4' 
on a workstation with Intel Xeon X5650 GPUs and one 
NVIDIA Tesla G2075 GPU card (late 2011). The second 
case considers the antibody immunoglobulin G, which is 
a much larger molecule than protein GB1D4'. For these 
runs, we used either: (1) Boston University’s bungee 
cluster, which has 16 nodes with 8 Intel Xeon GPU cores 
each, and a total of 3 NVIDIA Tesla K20 (Kepler, late 
2012) and 26 NVIDIA Tesla M2070/2075 GPUs; or (2) 
the George Washington University’s Colonial One clus- 


A. First case: protein G B1 D4' 


We investigated the preferred orientation of protein 
GB1D4' placed 2A away from a lOOAx lOOAx lOA block 
with surface charge density ±0.05C/m^, centered with 
respect to a lOOAxlOOA face. In biosensors, protein 
GBl D4' can be used as an intermediate protein, coupled 
to the functionalized surface directly by covalent bond¬ 
ing. The protein will thus be at a small distance from 
the surface. In this case, 2A is in the order of magnitude 
of the size of a water molecule, or of a C-N bond. The 
charge density of ±0.05G/m^ matches that used in other 
works .1^ In these cases, we considered a solvent with no 
salt, i.e., K = 0 (to compare with other published results), 
and with relative permittivity 80. The region inside the 
protein had a relative permittivity of 4. 

As seen in Figure atiit is the angle between the pro¬ 
tein’s dipole moment and the normal vector to the sur¬ 
face, and Q^rot rotates about the dipole moment. When 
the dipole-moment vector and the normal are aligned 
(u^tiit = 0)5 we define a vector Vref as the shortest dis¬ 
tance between the axis normal to the surface that goes 
through the center of mass, and the atom that is furthest 
away from it. We use Vref as a reference to define the ro¬ 
tation angle arot- the angle between Vref and the vector 
normal to a lOOAxlOA face. 

We sampled the total free energy every Ao^tiit = 2° of 
tilt angle and Ao^rot = 10° of rotation angle, resulting in 
3, 240 independent runs. The surface mesh had 4 trian¬ 
gles per square Angstrom on the protein geometry and 2 
triangles per square Angstrom on the charged surface. 

Numerical parameters are presented in Table [T| In a 
companion publication,^!^ we present a grid-convergence 
study using both an analytical solution and a case with 
protein GB1D4'.I^ We computed an approximate ex¬ 
act value of —222.43[kcal/mol] for solvation energy and 
317.98[kcal/mol] for surface energy using Richardson ex¬ 
trapolation with very fine parameters. With results that 
are less than 2% away from the approximate exact val¬ 
ues, we are comfortable with the parameters in Table [T| 
and mesh densities of 4 elements per square Angstrom 
on the protein and 2 elements per square Angstrom on 
the surface 

Using total free energy as the input, the integrals of 
Equation (13) can be computed by means of the trape¬ 
zoidal rule. Eigure presents the probability of the 
protein orientation in terms of cos((atiit), in intervals of 
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TABLE I: Numerical parameters used for numerically 
probing the orientation of protein GB1D4'. 


# Gauss points: Treecode: GMRES: 

in-element close-by far-away Ncrit P 0 toL 

9 per side 19 1 300 4 0.5 1Q~^ 


A cos((atiit) = 0.005 (Fig. and Aati\t=2° (Fig. [^. 
Table |n| presents the average orientation < cos((atiity> 
for the surface having either positive or negative charge 
density, and Figure shows the electrostatic potential for 
the preferred orientation in each case. 

TABLE II: Average orientation. 


< cos(otiit) > 

Negative Positive 
-0.968 0.963 


B. Second case: immunoglobulin G 

We computed the electrostatic field of immunoglobulin 
G—a protein widely used in biosensors—interacting with 
a 250Ax250Ax10A block, varying the conditions of sur¬ 
face charge and salt concentration. The protein was cen¬ 
tered with respect to a 250Ax250A face, at a distance 5A 
above it. In fabrication, antibodies are usually immobi¬ 
lized on the biosensor surface via a cross-linker molecule, 
which we model here by increasing the distance from the 
surface. As before, the solvent had relative permittivity 
of 80 and the protein of 4. 

Grid-convergence study for immunoglobulin G — 
Since this was the first time we did calculations on im¬ 
munoglobulin G, we carried out a grid-convergence study 
to make sure the geometry was well resolved and to find 
adequate values of the simulation parameters for sam¬ 
pling different orientations. The error plotted in Figure 

is the relative difference between the energy obtained 
using PyGBe with each mesh density and the estimated 
exact value computed with Richardson extrapolation. 

In this case, we computed the solvation energy and 
surface energy of a system consisting of a surface with 
charge density 0.05C/m^ and a protein with Gtiit =31° 
and Grot = 130°. Using the results from runs with 
a mesh density of 2, 4, and 8 elements per square 
Angstrom, we added the solvation and surface energies, 
and used Richardson extrapolation to obtain a value of 
—2792.22[kcal/mol], and an observed order of conver¬ 
gence of 0.85. This is our reference to calculate the errors 
in Figure There is a slight deviation from the expected 
value of the observed order of convergence (1.0), which 
we attribute to the non-uniform mesh generated by msms. 
Even though the mesh density is on average doubled for 


each run, there is no guarantee that the refinement is ho¬ 
mogeneous throughout the whole molecular surface. The 


numerical parameters are presented in Table III 


TABLE III: Numerical parameters used in the 
grid-convergence study with immunoglobulin G. 


# Gauss points: Treecode: GMRES: 

in-element close-by far-away Ncrit P 0 tol. 

9 per side 19 1 1000 6 0.5 10~^ 


Probing orientation of immunoglobulin G — We sam¬ 
pled the total free energy every A^tiit = 4° of tilt angle 
and A^rot = 20° of rotation angle, resulting in a to¬ 
tal of 810 runs. The surface meshes had 2 triangles per 
square Angstrom throughout. Numerical parameters are 
presented in Table m 

TABLE IV: Numerical parameters used in the runs 
probing orientation of immunoglobulin G. 


# Gauss points: Treecode: GMRES: 

in-element close-by far-away Ncrit P 0 tol. 

9 per side 19 1 300 2 0.5 10~^ 


With the computed total free energy, we obtained the 
probability of each orientation using Equation and 
the trapezoidal rule. We sampled all combinations with 
surface charges of a = ±0.05C/m^ and a = ± O.lC/m^ 
and salt concentrations of 145mM {k = 0.125 A“^) and 
37mM {k = 0.0625 A“^). For each of these cases. Figures 
and show a color plot of the probability distribution 
with respect to the tilt and rotation angles, and a 3D plot 
of the preferred orientation, where the solvent-excluded 
surface is colored by the electrostatic potential. 


C. Reproducibility and data management 

We have a consistent reproducibility practice that in¬ 
cludes releasing code and data associated with a pub¬ 
lication. The PyGBe code was released at the time of 
submitting our previous publication,^^ under an MIT 
open-source license, and we maintain a version-control 
repository. As with our previous paper, we also release 
with this work all of the data needed to run the numeri¬ 
cal experiments reported here, including running scripts 
and post-processing code in Python for producing the 
figures. To support our open-science goals, we prepared 
such a ^^reproducibility package^^ for each of the results 
presented in Figures UlIZl and the probability plots in 
Figures and The included running scripts invoke 
the PyGBe code with the correct input data and meshes 
(also included), and post-process the results to give the 
final figure, all with just one command. Please see the 
respective captions for a reference to the reproducibility 
packages, hosted on the figshare repository. 
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FIG. 5: Orientation probability distribution of protein GB1D4'. Figures |5a| and 5b are the probability with respect 
to the tilt angle and its cosine, respectively. Figures [5^ and [5d] are the probability distributions with respect to both 
the tilt and rotation angles. Data sets, figure files and running/plotting scripts are available under CC-BY.I^ 



(a) Negative surface charge (atilt = 172°, arot = 110°) (b) Positive surface charge (atut = 8°, arot = 150°) 

FIG. 6: Electrostatic potential of protein GB1D4' for the preferred orientations according to Figure Black arrow 
indicates direction of dipole-moment vector. 


V. DISCUSSION 

A. First case: protein G B1 D4' 

The orientation of protein GB1D4' near charged sur¬ 
faces was studied using a combined Monte Garlo and 


molecular dynamics method by Liu and co- worker^ and 
experimentally by Baio and co-workers.l^The availability 
of these published results was a motivation to use this 
protein for a first test, to compare with the results ob¬ 
tained with our model. 

The results presented in Figure show that for the 































FIG. 7: Grid-convergence study of the solvation plus 
surface energy for immunoglobulin G interacting with a 
surface with charge density of 0.05G/m^. Data sets, 
figure files and plotting scripts are available under 
CC-BY.l^ 


most likely orientations, the dipole-moment vector is 
aligned with the vector normal to the interacting surface. 
This indicates that the dipole moment is the dominant 
effect that determines the protein’s orientation, over local 
protein-surface interactions. This is the expected result, 
since protein GB1D4' is a relatively small biomolecule. 

Moreover, Figure reveals that protein GB1D4' be¬ 
haves like a point dipole, as the most likely orientations 
shift 180° when the sign of the surface charge is flipped. 
This is also explained by the dipole moment dominat¬ 
ing the orientation. In fact, we repeated this whole set 
of calculations but placing protein GB1D4' at a greater 
distance, SA away from the surface, and the results did 
not vary. 

The dipolar behavior described by our results agrees 
with the experiments done by Baio and co-workers,l3 
in which they observed opposite orientations of protein 
GB1D4' adsorbed on NH^ and GOO“ self-assembled 
monolayers. With positively charged surfaces, most of 
the proteins oriented with the N-terminal of the pro¬ 
tein pointing away from the surface, while for nega¬ 
tively charged surfaces the opposite occurred, with the 
C-terminal pointing away from the surface. This agrees 
with our results in Figure (the dipole moment vector 
of protein GB1D4' points from the G-terminal to the N- 
terminal). 

Liu and co- worker^ used a combined Monte Garlo and 
molecular dynamics method to obtain < cos((atiit) >= 
0.95 for cr = 0.05C/m^, and < cos((atiit) >= —0.85T0.05 
for (7 = — 0.05C/m^, which agrees well with our re¬ 
sults in Table nn Note that MD simulations consider 
van der Waals interactions and conformational changes of 
the protein, whereas these are not considered in our ap¬ 
proach, explaining the slight differences in < cos((atiit) >• 
However, as noted by other researchers ,1^^^ electrostatic 
effects often dominate protein-surface interactions and 
drive orientation during adsorption, while van der Waals 


effects play a role only in cases of very low surface charge. 
For example, in Ref. [H van der Waals effects were of 
consequence in a setup with surface charge of 0.006G/m^ 
and high ionic strength, leading to weak electrostatics. 
In a biosensor-fabrication scenario, this would only be 
the case with low-quality SAMs. 

The results with protein GB1D4' mean that an elec¬ 
trostatic solver with implicit solvent using the Poisson- 
Boltzmann equation is capable of capturing the driving 
mechanism of physical adsorption and orientation of the 
adsorbed molecule, at least in cases where the molecule’s 
dipole moment is dominating the orientation. This is im¬ 
portant because protein adsorption, being a free energy- 
driven process, is difficult to study experimental!}® and 
thus simulations offer a promising alternative. Full 
atomistic molecular dynamics, however, demands large 
amounts of computing effort, and the possibility of using 
an electrostatics solver may extend the range of systems 
that can be investigated. 


B. Second case: immunoglobulin G 

With our numerical model already verified using an 
analytical solution for spherical geometr}® and the suc¬ 
cessful results for protein orientation of a small protein 
near a charged surface (Section |V A| ), we proceeded to 
study the effect of surface charge and salt concentra¬ 
tion on the orientation of the antibody immunoglobu¬ 
lin G. Antibodies are widely used in biosensors as lig¬ 
and molecules, due to their affinity and specificity with 
the target molecule (antigen), and it is vitally important 
that they are adsorbed on the sensor with the antigen¬ 
binding Ig fragment (Fab) pointing away from the sensor, 
into the oncoming flow containing the antigens (known 
as “end-on” or “tail-on” orientation). Early experimen¬ 
tal studies found that antigen/antibody ratio was espe¬ 
cially low on negatively charged surfaces,!^ leading to the 
notion that protein orientation was affected to leading 
order by charge. One subsequent stud}® investigated 
the orientations of two iso-types of immunoglobulin G— 
IgGl, corresponding to PDB structure IIGY, and IgG2a, 
corresponding to PDB lIGT—adsorbed on positive and 
negatively charged surfaces. As an indirect method of 
probing antibody orientation, the researchers obtained 
adsorbed amounts and antigen/antibody ratios by means 
of surface-plasmon resonance experiments (e.g., a higher 
antigen/antibody ratio would indicate that more active 
sites are accessible and more antibodies are in a favor¬ 
able orientation). The finding was that IgGl mainly had 
a “head-on” (unfavorable) orientation on the negatively 
charged surfaces and a mix of “tail-on” (most favorable) 
and “side-on” orientations on the positively charged sur¬ 
faces. IgG2a, on the other hand, had many orientations 
on both surfaces with positive and negative charge, lead¬ 
ing to the conclusion that IgG2a is harder to control using 
electrostatic effects. Results consistent with these were 
obtained by Zhou and co- worker^ using a united-residue 
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(a) Probability for a = —0.05C/m^, k = 0.125A ^ (b) x-y plane view for atilt = 116° and arot = 100° 





FIG. 8: Orientation probability distribution and surface potential of the preferred orientation for immunoglobulin G 
near a negative surface charge. The black arrow indicates the direction of the dipole moment, and the circles enclose 
the Fab fragments. Data sets, figure files and plotting scripts available under CC-BY.^SI 
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(c) Probability for (7=0.1C/m^ and /^=0.125A ^ (d) y-z plane view for atnt = 32° and arot = 100° 
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(e) Probability for c7=0.05C/m^ and 6:=0.0625A ^ (f) y-z plane view for atut = 44° and arot = 120° 
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FIG. 9: Orientation probability distribution and surface potential of the preferred orientation for immunoglobulin G 
near a positive surface charge. The black arrow indicates the direction of the dipole moment, and the circles enclose 
the Fab fragments. Data sets, figure files and plotting scripts available under CC-BY.I^ 
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model: a coarse-grained model where each amino-acid is 
treated as a sphere. They find that IgGl will have the fa¬ 
vorable “end-on” orientation on positive surfaces, as long 
as the charge density was large enough (0.018C/m^, in 
their case) and the ionic strength was low. But IgG2a did 
not show a clear preferred orientation at the conditions 
they looked at; the authors attribute this to the weaker 
dipole moment of this iso-type. 

We investigated the orientation of IgG2a, which other 
studies found harder to orient favorably on a biosen¬ 
sor surface, and used two values of the surface charge 
(a = 0.05 and O.lC/m^) and two values of salt concen¬ 
tration (k = 0.125 and 0.0625A“^), in each case varying 
two-fold. Figures and present the probability dis¬ 
tribution of IgG2a for many orientations (given by o^tiit 
and Q^rot)? ill each case. The following discussion refers 
to each variation of the parameters and the effect on the 
preferred orientation of the adsorbed antibody and its 
probability. 


Effect of surface charge — The lower value of sur¬ 
face charge here is a = ±0.05C/m^, the same value 
used in Ref. |8] to mimic the experiments reported in 
Ref. [7| Figures and show that for the lower 
value of surface charge with the higher salt concentra¬ 
tion [k = 0.125A“^), there is no clear preferred orienta¬ 
tion, to the point that the highest probability falls under 
10%. This means that adsorbing the antibodies under 
these conditions would result in a wide range of orienta¬ 
tions, which would not be favorable for biosensor fabri¬ 
cation. Moreover, the preferred configurations in figures 


8a and 9a show the antibody lying flat on the surface. 


far from the desired “tail on” orientation. This observa¬ 
tion is consistent with a previous study using a unified- 
residue model,!^ where this particular antibody showed 
many possible orientations. The authors of that study 
attributed this behavior to the weaker dipole moment of 
this molecule, compared with the variant IgGl. 

With the higher value of surface charge, in this case 
(j = ±0.1C/m^, the orientation probability distribution 
in the case of negative charge improves somewhat, as the 
antibody is slanted sideways rather than lying down for 
n = 0.125A“^ (at least one antingen-binding fragment 
is pointing up), and the probability of the preferred ori¬ 
entation is almost doubled for low salt concentration, in 
a “side on” orientation. For positive surface charge the 
slanted orientation is similar, however the probability is 
higher for the preferred orientation in both the low- and 
high-salt cases. In the cases with higher salt concen¬ 
tration, Figure shows a preferred orientation with a 
higher probability of 12%, compared to 8% in Figure 
and the dipole moment rotates towards the normal vec¬ 
tor. For the lower value of salt concentration (Figure [^, 
this effect is smaller, however it shifts the preferred tilt 
angle in the opposite direction, from 44° to 64°. Note 
that the dipole-moment vector does not point straight 
through the middle between the two Fab fragments, but 
in an angle. This indicates that, in contrast to protein 
GBl D4', local interactions dominate over the dipole mo¬ 


ment. If the dipole moment were the dominant effect, 
the dipole-moment vector would tend to align to the sur¬ 
face normal as the surface charge increases. Th is argues 
against the suggestion by other researcher^^^^ that the 
dipole-moment vector is the main determinant of orien¬ 
tation. 


Effect of salt concentration — As the surface charge 
density was varied two-fold, we also varied the Debye 
length (/^“^) two-fold. In terms of salt concentration, it 
means a 4x decrease in the amount of salt. The higher 
value of salt concentration corresponds to 145mM, which 
is in the physiological salt range. 

Like increasing the surface charge, lowering the salt 
concentration affects the orientation probability distri¬ 
bution. For (j = —0.05C/m^ (Fig. [^, the effect is a 
large shift in the preferred tilt angle, from amt = 116° to 
= 40°, with a small change in probability. For the 


positive weaker charge, a = 0.05C/m^ (Fig. 9e), not only 
does the peak probability increase considerably (^ 3x), 
but the preferred tilt shifts from 64° to 44°. This ori¬ 
entation is favorable for biosensing applications, as the 
antigen-binding fragments are pointing away from the 
surface, in a “tail on” orientation. For the stronger neg¬ 
ative charge, a = — O.lC/m^, the probability peak in¬ 
creases 2.5X for a “side on” orientation where one of the 
antigen-binding fragments is attached to the surface in 
an unfavorable position (Fig. [^. With positive surface 
charge, the tilt angle shifts in such a way that the anti¬ 
body is lying on the surface with a marked probability 
close to 30% (Fig. 9g). This orientation is not ideal for 
biosensors, but it is better than the slanted position as 
neither of the Fabs are attached to the surface. 

From the results in Figures and we conclude that 
the iso-type IgG2a can in general be better orientated 
with low salt concentration and high surface charge, as we 
get more pronounced high-probability regions. Moreover, 
good orientations for biosensors are more likely to occur 
with positive surface charge (Figure [M|, since the Fab 
fragments are pointing up. Previous studies had shown 
that the IgGl variant could be controlled, but not IgG2a. 
The advantage of a positive surface charge and a low ionic 
strength had been suggested by previous studies, but not 
for this particular variant of immunoglobulin G. Note also 
that our lower value of salt concentration is 3 7mM, which 
is a higher amount of salt than other studies 


A limitation of this study stems from the appli¬ 
cation of linearized Poisson-Boltzmann equation. Bn 
and co-workerPS assessed the accuracy of the Poisson- 
Boltzmann equation for highly charged surfaces (^ 
0.4C/m^), getting good agreement of the model with ex¬ 
periments. Rigorously, the linearized Poisson-Boltzmann 
equation is a valid approximation of the Poisson- 
Boltzmann equation when the nondimensional potential 
is smaller than 1 {f^qe/ksT << 1). However, this restric¬ 
tion can be relaxed when calculating solvation energy. 
For example, we ran a calculation using our boundary 
element code on an isolated IgG2a immersed in a solvent 
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with 37mM of salt {k = 0.0625A“^), with the parameters 
from Table [IV| Computing the absolute value of the di¬ 
mensionless potential on the molecular surface gives over 
55% of the triangles with (^qe/ksT > 1 and an average 
value of 1.5. Yet, using the linear or non-linear Poisson- 
Boltzmann model of APBS,^ the solvation energy for 
the isolated IgG2a gives the same result. (The APBS 
tests used a volumetric mesh of 150 x 150 x 150A^ with 
449^ nodes.) Adding a surface with charge cr = O.lC/m^ 
next to IgG2a in the configuration of Fig. the situa¬ 
tion is similar: 51% of the triangles have a dimensionless 
potential with absolute value exceeding 1 and the aver¬ 
age is 1.3. By comparison with the isolated IgG2 case, 
we expect linearized Poisson Boltzmann to give a good 
approximation of the solvation energy in this case. This 
conclusion is also consistent with the results of Ref. |46| 
that show good agreement between linear and non-linear 
Poisson Boltzmann when the average dimensionless po¬ 
tential on the molecular surface is between —2 and 2. 
Finally, with a value of surface charge of O.lC/m^, the 
mechanism of protein orientation appears to be dictated 
by the solvation energy, rather than the surface energy. 
That is, the maximum probability of the preferred orien¬ 
tation occurs at the minimum of solvation energy.l^ Since 
solvation energy in our systems is well approximated by 
linear PB, we conclude that the use of this theory is jus¬ 
tified in these cases. 

Our results suggest that a combination of high surface 
charge and low salt concentration increases the proba¬ 
bility of the preferred configuration, and that more fa¬ 
vorable orientations are obtained with positive surface 
charge. We completed an additional set of tests for the 
orientation of IgG2a near a surface with a = 0.2C/m^ 
and K. = 0.03125A“^, even though the electrostatic po¬ 
tential is outside the linear regime in this case. The result 
in Figure [^suggests that as we increase the electrostatic 
effects, the preferred orientation becomes highly marked, 
and tends towards a favorable orientation for biosensors. 
However, our model is based on the linearized Poisson- 
Boltzmann equation and this test goes outside the known 
regime where the model can be applied. We cannot claim 
that this is a physical situation, but it indicates a trend. 
Further study of the conditions that make IgG2a orient 
favorably for biosensors may need nonlinear models, and 
combined experimental observations. 


VI. CONCLUSION 

Various studies have revealed the importance of pro¬ 
tein orientation in immunoassays. One work suggested 
that highly oriented antibodies could result in 100 x im¬ 
provement in the affinity of a biosensor .1^ Thus, a design 
goal would be to know how to prepare a surface to control 
protein orientation. Yet, despite much work, control of 
protein orientation has not been successful. This study 
increases our understanding of how nanosurface proper¬ 
ties (charge) and preparation conditions (salt levels) af¬ 


fect protein orientation. We successfully used an implicit- 
solvent model to study protein orientation near charged 
surfaces, which in our method can have any geometry. 
In a companion publication,^ we describe expanding the 
applicability of our open-source code, PyGBe, to account 
for the presence of charged surfaces and present grid- 
convergence studies using an analytical solution and pro¬ 
tein GB1D4'. 

Protein GB1D4' behaves like a point dipole near a 
charged surface, with the dipole-moment vector shifting 
^180° when the sign of the surface charge flips. Our re¬ 
sults compare well with experimental observations and 
simulations using combined Monte Carlo and molecu¬ 
lar dynamics methods, supporting the use of our ap¬ 
proach for probing protein orientation near charged sur¬ 
faces. We applied our approach to immunoglobulin G, 
a biomolecule that is much larger than protein GB1D4' 
(about 125 X, by volume) and would be challenging to 
study via molecular dynamics. The iso-type IgG2a was 
found by previous studies to be hard to control, exhibit¬ 
ing many orientations, but we are able to obtain a pre¬ 
ferred orientation that is favorable for biosensing with a 
positive surface of 0.05C/m^ or higher d 37mM of salt 
in the solvent. We conclude that local electrostatic in¬ 
teractions dominate over the dipole moment, and even 
this protein can be favorably oriented with the appropri¬ 
ate fabrication protocol. Potentially, protein engineering 
could be used to obtain ligand molecules that interact 
with charged surfaces in a desired fashion. In this appli¬ 
cation, where ligand molecules undergo little conforma¬ 
tional change as they adsorb on the sensor surface, our 
new implicit-solvent model can offer a valuable approach 
to assist in biosensor design. In our future work, and in 
collaboration with experimental researchers, we intend 
to use this approach to aid the design of better ligand 
molecules, by looking at the preferred orientations for 
different ligand mutants. 
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